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ABSTRACT 

The Schrodinger equation was solved numerically for a one- dimensional collision be- 
tween a harmonic oscillator and a particle, with a rectangular trapping potential. Gener- 
al solutions were sums of states, each a product of a particle wave function and an oscil- 
lator wave function. A correlation was suggested between probability current densities 
for states within the well and sticking of gas atoms on surfaces, where sticking is regard- 
ed as arising from multiple collisions. Energy accommodation coefficients were formu- 
lated and compared with helium- tungsten experiments for the temperature range 20 to 
120 K. The theory underestimates the accommodation coefficient at lower temperatures 
and overestimates it at higher temperatures. Results indicate that only approximate 
qualitative predictions of the coefficient are possible with this simplified model. 
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QUANTUM THEORY OF A PARTICLE-OSCILLATOR COLLISION APPLIED 
TO GAS-SURFACE INTERACTIONS AT CRYOGENIC TEMPERATURES 

by Jerry D. Smith 
Lewis Research Center 

SUMMARY 

A mathematical model of the interaction of an inert gas atom with a clean cold sur- 
face of a solid is proposed. The model is a one-dimensional collision between a particle 
acting as the gas atom and a simple harmonic oscillator acting as a surface atom. The 
interaction potential between the two is taken as hard sphere, with a rectangular well in 
front of the infinite wall. Two regions are created: the well region, where the particle- 
oscillator separation distance is less than a certain value, and the region outside the well, 
where the separation exceeds that value. 

The solutions to the Schrodinger equations in the two regions are written as sums of 
component states. Each such state is a product of an unknown coefficient, a particle ex- 
ponential wave function, and a harmonic oscillator wave function. Boundary conditions 
are applied, and the unknown coefficients are derived numerically. Probability current 
densities for component states neighboring a low initial oscillator state are calculated in 
the two regions for helium and neon on tungsten. 

A gas-surface interaction is simulated by an ensemble of particle-oscillator systems. 
A correlation is suggested between the sum of probability current densities within the well 
and sticking of gas atoms on a surface. This sum is physical mass flux in one direction. 
We interpret sticking as arising from a gas atom making multiple collisions with a surface 
atom. An increase of particle flux within the wells is indicated by the probability current 
density sum within the wells exceeding the initial incoming density. 

Outside the well, transition probabilities are derived from current densities. Then, 
microscopic energy accommodation coefficients are formulated, suitably averaged, and 
compared with helium- tungsten experiments for the temperature range 20 to 120 K. The 
theory underestimates the accommodation coefficient at the lower temperatures and over- 
estimates it at the higher temperatures. Results indicate that only approximate qualita- 
tive predictions of the coefficient are possible with this simplified model. 


INTRODUCTION 


Variations of the present particle-oscillator collision model have been applied before 
to rarefied gas-surface interactions. The particle-oscillator collision is a reasonable be- 
ginning because, for free molecular flow, momentum and energy transfer in a gas-surface 
interaction depend entirely on the nature of the collisions of individual gas atoms with sur- 
face atoms (ref. 1). The model in the present treatment is a one-dimensional collision be- 
tween a simple harmonic oscillator simulating one of the atoms comprising the surface of 
the solid, and a particle acting as the gas atom, with a hard sphere interaction potential 
modified by a rectangular well between the oscillator and the particle. 

This model and its mathematical treatment are based on the treatment given by Shuler 
and Zwanzig (ref. 2) for a similar model. They used only a hard sphere interaction poten- 
tial (no well) and calculated oscillator transition probabilities. Using a computer solution 
of the Schrodinger equation, they showed that multiphonon energy exchange can have ap- 
preciable probability for certain initial incoming particle energies. This contradicts the 
one-phonon-selection rule evolved in perturbation treatments of particle- oscillator models 
(e. g. , Devonshire, ref. 3). This one-phonon rule has been shown inadequate to explain 
condensation of gases on surfaces (ref. 2). 

Gilbey used the same model as Shuler and Zwanzig, in a classical and quantum treat- 
ment of gas-solid energy exchange (ref. 4). His calculation procedure necessitated a 
choice of low particle-to-oscillator mass ratio and high particle-to-oscillator energy ra- 
tio. Sec rest and Johnson used an exponential repulsion interaction potential in their col- 
lision model (ref. 5). They calculated transition probabilities using the method of ampli- 
tude density functions. Jackson and Mott (ref. 6), Castellan and Hulburt (ref. 7), and 
Widom (ref. 8) give special approximate analytical solutions for impulsive inelastic colli- 
sions. 

Of the previous models, only Devonshire's (ref. 3) uses a trapping interaction poten- 
tial which allows the particle to be found in a bound state near the oscillator. Devonshire's 
model uses a Morse potential. The previous treatments, including Devonshire's, yield 
only asymptotic forms of the wave function as the particle approaches +°® after the inter- 
action. The objective of these previous treatments was to formulate transition probabil- 
ities or accommodation coefficients; therefore, primary interest lay in the final state of 
the system. In contrast, the present treatment yields solutions in a well region near the 
oscillator. And it does not require assumptions in the mathematical development either 
about the relative magnitudes of the particle and oscillator coordinates, or about relative 
magnitudes of initial particle energy, initial oscillator energy, and interaction potential 
well depth. 

Diestler and McKoy (ref. 9) recently presented a general method for the quantum- 
mechanical treatment of the inelastic collision of composite particles. Their method 
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could include the Shuler- Zwanzig (ref. 2) and Secrest- Johnson (ref. 5) models, among 
others, as special cases. However, to adapt the formalism of Diestler and McKoy to the 
present model, one would have to require that (1) the total wave function vanish when the 
oscillator coordinate exceeds some particular value, and (2) that the interaction potential 
vanish when the particle coordinate exceeds some other particular value. Condition (2) 
defines their asymptotic region. The interaction potential of the present report vanishes 
when the difference between the oscillator and particle coordinates exceeds some given 
value. And condition (1) is not required in the present treatment. Furthermore, and more 
important, although Diestler and McKoy can use a trapping interaction potential, they can 
derive Schrodinger equation solutions only in their asymptotic region, not in our interac- 
tion region (the rectangular well). 

Sugawara gave a quantum- mechanical theory of adsorption using basically the same 
model proposed in the present study (ref. 10). But he allowed a single group of phonons 
of identical energy to leave the particle-oscillator system at a definite time during the in- 
teraction. He then derives an expression for the probability of finding the gas particle in 
a trapped state immediately after this energy release. However, he was not able to veri- 
fy, with any experimental comparisons, that either this mechanism of energy release to 
the crystal or the various approximations introduced in his time-dependent perturbation 
treatment were justified. 

The present model also does not now yield sticking coefficients to be compared with 
experiments. We only propose a different outlook on the quantum mechanics of trapping 
a gas atom on a surface: linking the multiple collisions that occur when a gas particle re- 
sides in the vicinity of a surface atom to the particle probability current density in that 
vicinity. The model, however, does yield simply formulated energy accommodation coef- 
ficients. But the results for helium on tungsten indicate that only approximate qualitative 
predictions can be made with this simplified model. 


MODEL 

The incident gas atom interacts along the line of centers with a rotationless harmonic 
oscillator, as shown in figure 1. There are several potentials which could be used to de- 
scribe the physical interaction between the particle and the oscillator. Such potentials 
are attractive when the interatomic separation is large and repulsive when the separation 
is small. A hard-sphere potential modified by a rectangular well in front of the infinite 
wall (see fig. 2) was chosen because then the Schrodinger equation remains separable in 
the interaction region. 

In appendix B, values of the well width b and well depth V ^ were chosen so that the 
rectangular well potential would approximate a Morse potential. And for this approxima- 
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Figure 1. - Particle-oscillator interaction. Figured - Rectangular interaction potential V Mm . Coordinate r is 

linearly related to particle coordinate r 1 (see appendix B), and R is 
oscillator coordinate. 


tion, the Morse potential parameters given by Goodman for rare gas atom - metallic atom 
combinations were used (ref. 11). The coordinate r in figure 2 is linearly related to the 
actual particle coordinate r' and is introduced to facilitate subsequent equations (see ap- 
pendix B). The Schrodinger equations are invariant to this linear transformation. 

In terms of the coordinates r and R, the Schrodinger equation in region 1 (r - R 
> b) is 


Mco 2 R 2 - iL-J^W (r,R) = E(l)*.(r,R) 
\ 2M 3R 2 2 2m 9r 2 / 


( 1 ) 


and in region 2 (0 ^ r - R < b) is 


2 2 2 2 

2L J_ + 1 Mgo 2 R 2 - — - V. l* 9 (r, R) = E(2)tf 9 (r, R) 

2M .„2 2 2m . 2 1,2 1 


( 2 ) 


9R 


9r 


All symbols are defined in appendix A. 

Equations (1) and (2) are separable. The most general solutions to these equations 
are linear combinations of harmonic oscillator wave functions, ?// ( R ), with oscillator en- 
ergy E qsc = ^n + -^hu>, n = 0, 1, . . . , and free particle wave functions: 


^(r.R) = A^j(R)exp|- jfepj + Y Bji//j(R)expi+ 

f=0 ^ 


i[Pl(f)] 


r) 


K 


(3) 
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( 4 ) 



The indices j and f denote initial and final states of the oscillator, respectively. Parti- 
cle momentum is given by energy conservation in the two regions: 


Region 1: 


Region 2: 



(5) 


( 6 ) 


Now, in the total wave functions, states for which Pj(f) and p 2 (f) are imaginary are 
included. These nonphysical states arise when f is large, as can be seen from equa- 
tions (5) and (6). They must be included in the most general solutions, because without 
them all the boundary conditions cannot be satisfied. Briefly, the omission of such states 
would lead to an infinite number of linear equations in a finite number of unknowns, which 
cannot, in general, be solved. (This set of equations would differ from eqs. (17) to (20) 
only by having finite summations over states f, instead of infinite summations. ) 

The notation Pj(f) and p 2 (f) is used to indicate that, in calculating p 1 (f) and p 2 (f) 
from equations (5) and (6), the positive root is always chosen, whether it is real or imag- 
inary. Doing this ensures that each possible component state is included in the general 
solutions just once. The initial state appears as the first term in ^(t, R). Writing the 
general solutions as shown herein makes the physical meaning of the various component 
states clear; particle motion to the left or right is implied by negative and positive values 
of real momentum, respectively. 

The boundary conditions are 


* 2 (R,R) = 0 


(7) 


5 


r-R=b 


( 8 ) 


— *i(r,R) 
0R 1 


r-R=b 


= J-* 2 (r,R) 
3R 1 


d_ 

0r 


^(r, R) 


r-R=b 


— * 2 ( r > R ) 
9r z 


r-R=b 


(9) 


* x (r, R) 


r-R=b 


= * 2 (r, R > 


r-R=b 


( 10 ) 


The first condition (eq. (7)) is the only one that can be imposed properly at the infinite 
wall. 

It is important to note that the boundary conditions are imposed on the total system 
wave functions. The only energy quantization that occurs is the well-known E Qgc = 

^n + as oscillator energy. This quantization arises only be requiring the wave func- 

tions to be bounded everywhere and approach zero as R — °°. The only requirement for 
particle momentum is that the corresponding energy satisfy energy conservation. To 
evolve the familiar energy levels for a particle in a well, one would have to impose func- 
tion and derivative continuity on parts of the total system wave functions, namely, the 
free particle exponentials, and this is incorrect. Since the total system Hamiltonians in 
both regions do not depend on time explicitly, a time-independent formulation is possible, 
and the total system energy is constant. However, because of energy exchanges with the 
oscillator, particle energy is not constant in time. Therefore, previous elementary time- 
independent particle- in-w ell treatments are inapplicable. 

The mathematical procedure introduced by Jackson and Mott (ref. 6) and used by 
Shuler and Zwanzig (ref. 2) will be used here. The general objective is to evolve a set of 
simultaneous linear equations which can be solved by computer for accurate values of the 
lower ordered unknowns B^, C^, and Dj (i. e. , those having low values of f). These val- 
ues are used to calculate probability current densities for the lower ordered states: 


^| Cf | 2 


P2< f > 

m 



2 


6 


The oscillator will be given an initial low state, j = 0, 1, 2, and sometimes 3, approxi- 
mately corresponding to a low wall temperature, as seen from appendix E. Interest 
naturally centers on transitions to neighboring low states. The initial incoming density 
is [p(j)/m]| A.| 2 , and A- is taken equal to 1. 

It will be convenient to nondimensionalize by expressing all energies in terms of fiw 
and all lengths in terms of the characteristic harmonic-oscillator length a~^ = (Mw/E)" 
by using the following relations: 



Kco 


*\ 


fi = ba 


X = Ra 


2 2m 

r = — 

M 


G 0 " 


Ml! 

2mKco 


£ i = 




2mKco 


= e n + j - f 


e 2 " 


[ P 2<«] 2 _ 


2mKco 


- Cq + j - f + fi 


- y £ 


1/2 

0 


( 11 ) 


h l = r< c 0 + J ■ f ) 1//2 


y(e 0 + i - f + /3) 


1/2 


As previously noted, because of the form of the general solution, the positive root is al- 
ways chosen in calculating h. or h 0 . For example, if for some f, e. < 0, h. is taken 
to be ivlcjl 1 / 2 . 12 11 

Using equation (4) and nondimensionalizing, equation (7) becomes 


0 j^(X)|Cj exp(-ih2X) + Dj exp^ihgX)! (12) 

f=0 


Some terms in equation (12) can be expanded in the complete orthonormal set of harmonic 
oscillator wave functions ji// n (X)J: 


exp(-ih 2 X)ify f (X) P K (h 2 ,f W K (X) (13) 

K=0 
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( 14 ) 


OO 


exp (-nhgX^X) =£ Q K (h 2 ,f)^ K (X) 
K=0 


The matrix elements involved are readily available (ref. 7): 


Q K (h2, f) - 



1/2 




(-iMhg )- 21 
t:(K - t)i(f - 1): 


(15) 


where t is the smaller of the values f or K. When +i is replaced by -i, P„(h 0 ,f) 
in z 

is identical to Q K (h 2 , f). Substituting these expansions into equation (12) yields 


0=E ^K (X )| 

K=0 


E C f P K (h 2 ,f)+D f Q K (h 2 ,f) 


f=0 


(16) 


Now, since |^(X) j is complete and orthonormal, every coefficient of i// K (X) in equa- 
tion (16) must vanish: 


0 - E [ C f P K< h 2- £ > + WV £ >] 

f-0 


K - 0, 1, ... , 


(17) 


As shown in appendix C, the other three boundary conditions (eqs. (8) to (10)) can be 
treated similarly, yielding the following three sets of equations: 

2|^j 6 fj exp(-ih oM )[-s K (h 0 ,j) + ( 2 i> 1/2p K (Vi ' X) ] 

+ B f exp(ih^ja)^-T^(h^,f) + (2f) l/2 Q K (h 1 ,f - 1)J 
+ C f exp(-ih 2M )[s K (h 2 ,f) - (2f) 1 / 2 P K (h 2 ,f - 1)] 

+ D f exp(ih 2fl ) [T K (h 2 , f) - (2f) l/2 Q K (h 2 , f - 1)J | = 0 K = 0, 1, . . . , (18) 
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K = 0, 1, ... 00 


(19) 


eo p 

2 L 6 fj A j h 0 e3 ^(- ih oM)P K (h 0 ,f) + B f hj exp(ih 1 n)Q K (h 1 ,f) 


f=0 


+ C f h 2 exp(-ih 2 /i)P K (h 2 ,£) - D f h 2 exp(ih 2 /i)Q K (h 2 , f)J = 0 
£ [ 6 fj A j exp(-ih 0 /i) p K ( h 0 , j) + B f exp(ih 1 p)Q K (h 1 , f) 


- C f exp(-ih 2 p)P K (h 2 , f) - D f exp(ih 2 p)Q K (h 2 , f)J = 0 K = 0, 1, ... 00 (20) 

CALCULATION PROCEDURE 
Numerical Approximation Technique 

This fourfold infinite number of equations (eqs. (17) to (20)) in a threefold infinite 

number of unknowns (B^, C^, D^, f = 0, 1, °°) is approximated by a finite order matrix 

equation, which is then solved by computer. The f summation is taken only to some in- 
teger F, and K is cut off at some integer K. Somewhat arbitrarily, each boundary con- 
dition is allowed to contribute equally to the approximating finite set by requiring that 

3(F + 1) = 4(K + 1) (21) 

The left side of equation (21) gives the total number of unknowns in the finite set, while 
the right side gives the total number of equations in that set. Taking F and K high 
enough gives sufficiently converged values of | cj and | D^| for small values 

of f. The computer program was written so that any positive integers F and K satis- 
fying equation (21) could be used, thereby enabling the convergence of the answers of in- 
terest to be easily tested. 

The input for the program is 

€q initial incoming particle energy, units of fico 
p nondimensional rectangular well width 

0 nondimensional rectangular well depth 

y mass ratio 

j initial oscillator state 

F,K cutoff indices 
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The output of the program included values of all the unknown coefficients Bq, Bp . . . , 
B f , Cq, Cp . , C F , Dq, Dp..., Dp separated into real and imaginary parts. 

Current densities for the various component states were normalized by dividing them 
by the initial incoming density [p(j)/m]| A . | with Aj = 1. These normalized densities 
become (in terms of dimensionless units) 



+ 3 - f + P 



2 


+ j-f + A l/2 , 


I D r 


re 0-- f \ 1/2 , Bf |2 

: 0 


and are referred to as C^- state densities, D^- state densities, and B^- state densities, re- 
spectively. For the cold-rare-gas-atom - cold- metallic-atom combinations of interest, 
eg, j, and /3 are all low, < 3 or 4. This leads to real, physically meaningful current 
densities only for low values of f, for which the squares of the coefficients | 2 , etc.j 
have accurately converged. When these densities were real, they were computed. Then, 
the normalized real C f -state densities were summed, and this sum is defined as a C- state 
sum. Similarly, the appropriate D^- state densities and B^- state densities were summed. 
Convergence of these three sums was ensured by convergence of the individual densities. 

Two other sums were formed by adding all normalized C^-state densities, and then 
Dj- state densities, for which particle energy is less than zero but above the well bottom. 
These sums are referred to as bound- state sums. 


Energy Accommodation Coefficient Calculation Procedure 


The energy accommodation coefficient is defined as 


a 


ac 


lim 

E r E o V E ° " 



( 22 ) 
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where Eq is the mean energy of the incident particle, is the surface molecular ener- 
gy, and E2 is the average energy of particles leaving the surface. The definition is often 
stated in terms of corresponding temperatures, and is strictly accurate only if the tem- 
perature Tg of particles leaving the surface is well defined. This is true if such parti- 
cles have a Maxwellian velocity distribution. This distribution will be undistorted only 
when Tq ~ Tj (ref. 12). The better experiments are conducted so that AT = Tq - T^ (or 
AE = Eq - Ej) is as small as possible. 

The computer results for the microscopic particle- oscillator interaction in the pres- 
ent model can be used to formulate theoretical values of a . A microscopic energy ac- 
commodation coefficient is defined as 


E. - E . 
a = in out 

3 E. - E 
in osc 


(23) 


where 


E osc 

E out 


incoming particle energy, CgRcj 
energy of oscillator, (] + 


average particle energy after interaction, calculated by weighting each of the pos- 
sible stationary-state particle energies in the asymptotic region by calculated 
transition probabilities to these states, [(e 0 + 3 “ f )A 0 )] 1/2 | B f! 2 = P j-f’ and 
then summing these weighted energies 


This Qfj is similar to Goodman's "effective a. c. " of a single gas atom (ref. 11). 

The surface is considered to be composed of independent oscillators, with the popu- 
lation of states given by the Boltzmann law 


P(j) = Fraction of oscillators in state j 
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where T is the effective surface temperature. If the incoming particle strikes the sur- 
s 

face head on at a random oscillator, P(j) is also the probability that the particle strikes a 
j-state oscillator. Each oscillator is assumed to occupy the same fraction of surface 
area. This surface model is admittedly a poor one. For example, particles that miss the 
oscillators cannot be accounted for. 

A thermally averaged accommodation coefficient « calc is defined by the equation 


“calc =E < 25 > 

i=o 

In this equation, P(j) is regarded as constant in time by regarding the surface tempera- 
ture T as being continually maintained by an external reservoir. A sample helium- 
tungsten calculation is given in appendix E, and the unimportance of states higher than 
j = 3 in equation (25) is shown. 


Accuracy Checks 


Several checks were made on the accuracy of the calculations. For the model of ref- 
erence 2 (identical to the present model with zero well depth), Shuler and Zwanzig calcu- 
lated probabilities of transitions in which the oscillator goes from state j to state f, 
while the incident particle momentum goes from -p^ to +p^. These probabilities, as in 
the present model, are given by 


' e n + 5 - f 


A/2 


|B f | =P M 


Their answers were duplicated when the well depth in the present model was set equal to 
zero, and y = 1. 00, for a range of well widths 0 < ju < 100. Computer results also indi- 
cate that with a depthless well, ^(r, R) = ^(r, R), as far as the accurately calculated 
coefficients B f , C f , and D f are concerned, as expected. 

Microscopic reversibility (a result from energy conservation) holds in the present 
model in the form 



e=e 


0 


= P 


f-J 


e=e 


o +j - 


(26) 
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Conservation of current density provides checks on the B-, C-, and D- state sums de- 
fined previously. The total wave functions (eqs. (3) and (4)) can be used to show that the 
particle conservation requirement (total particle current density at any point must be zero) 
implies that the B- state sum must be equal to 1, and that the C- and D- state sums must be 
equal to each other. 

The conventional form for current density is 

S = JL. (**v* - V***) = — Im(**V*) (27) 

2im m 

If particle motion without regard to the oscillator motion is of principal interest, the oscil- 
lator coordinate R in equation (27) is integrated out (ref. 13, pp. 126-128). Separability 
in the problem and normalization of the total wave functions with respect to R, but not r, 
make this the most straightforward procedure to follow to derive a relation satisfied by 
the coefficients. Setting the net particle current densities in the two regions equal to 
zero, one obtains 

/ OO 

o Im[^*(r,R)V r ^ 1 (r,R)]dR = 0 (28) 

y "0O 

Im[tf*(r, R)V * 2 (r, R)ldR = 0 (29) 

As shown in appendix D, equation (28) implies that the B-state sum must be equal to 
1, and equation (29) implies only that the C- and D-state sums must be equal to each 
other. The common value of the C- and D-state sums may be less than, equal to, or 
greater than 1, depending on the input parameters. (The variance of these sums from 
1 can be demonstrated analytically with a simplified version of the present model, having 
a potential like that in fig. 2 but fixed in space (e. g. , R is set s o) with no energy ex- 
change possible. ) A physical explanation for this difference between incoming particle 
current density and that within the well is given in the RESULTS AND DISCUSSION sec- 
tion. 


RESULTS AND DISCUSSION 

Proposed Application to Trapping of a Gas Atom on the Surface of a Solid 

Sticking probabilities are usually defined in experiments as fractions of incident gas 
atoms which are adsorbed on a surface, often by an expression of the form 
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Incident beam intensity - Outcoming beam intensity 
Incident beam intensity 

In many sticking probability experiments, the incident gas is not a beam, but just the part 
of a gas, injected into an evacuated chamber at a known rate, which strikes a known ad- 
sorbing surface. After the chamber pressure stabilizes, the adsorption rate on the sur- 
face must be just the flow rate into the chamber. Then, the sticking probability is just 

Rate of adsorption on surf ace 
Rate of collision with surface 

where the denominator is calculated from kinetic theory. 

In the present model, the total particle current density leaving the well equals the in- 
coming density (the B-state sum equals 1). Thus, it is difficult to discuss trapping by 
using this definition of sticking probability. 

However, in some classical models (e. g. , Goodman’s in ref. 11), the sticking proba- 
bility is taken as the fraction of gas particles which undergoes more than one collision 
with the surface, specifically, with a single surface layer atom. Recently (ref. 14), 
Goodman referred to such trapping as "apparent” trapping and stated that classical the- 
ories that use such a definition greatly overestimate the few reliable experimental meas- 
urements of trapping probabilities. However, he stated that "any one definition is prob- 
ably not the most suitable for all purposes. " These classical models allow energy to be 
exchanged between the oscillator-particle system and a one- or three-dimensional lattice 
during the interaction. This is not done in the present model. (Only in the energy accom- 
modation coefficient calculation were oscillators permitted to be restored to their original 
states, and then only after the interaction. ) 

In a quantum- mechanical model, conditions cannot be determined for a gas particle 
definitely to undergo more than one collision with the surface. Nor can meaning be as- 
cribed to a collision with the oscillator in a prescribed phase, or to an n th collision dur- 
ing the interaction. Although these restrictions are recognized, trapping in the present 
model can still be discussed from the standpoint of multiple collisions. It is suggested 
that probability current densities in the well may be interpreted physically in terms of 
multiple collisions. The conditions under which a gas particle is most likely to undergo 
more than one collision with the surface are determined, as shown when the results for 
helium and neon on tungsten are discussed. 

Current densities calculated within the well are shown by Beder (ref. 15) to be close 
numerical approximations to actual physical mass fluxes in the various states, as would 
be observed in an ensemble of classical systems. He gives an involved analysis of mass 
fluxes in a potential similar to that of figure 2, but fixed in space (not fastened to a moving 
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oscillator). Briefly, he states that if the particle energy eigenfunction is approximately 
an expansion in particle momentum eigenfunctions for a limited region of space, this ap- 
proximate expansion will furnish approximate relative particle number densities for the 
accessible momentum states, when measurements are taken anywhere in the limited re- 
gion. In the present model, the particle energy eigenfunction is immediately an expan- 
sion in particle momentum eigenfunctions (plane waves) and hence will furnish meaning- 
ful particle densities. 


Results for Helium and Neon on Tungsten 


The input parameters for helium on tungsten are as follows: 


y = 



0. 2087 


k$D 

= 0. 3 (using appendix B of ref. 11) 

K 


co = 1. 31X10 13 sec -1 


for tungsten 


0 = 0. 7045') 


M = 20. 76 


f(see appendix B) 


li = 36. 33 J 


The initial particle energy, in units of Kw, was varied from £q = 0. 025 to 1. 50 (with 
another calculation at €q = 1. 90). For the approximate correlation ^ kT — CgRco, this 
6 q range corresponds to a gas temperature range of 5 to 300 K. (The mean kinetic en- 
ergy of a one- dimensional Maxwellian beam is ^ kT (ref. 2, p. 1783). ) Temperature 
must be correlated with energy only for a one- dimensional beam to obey the principal of 
detailed balancing (ref. 2, p. 1783). This principal ensures thermal equilibrium in a 
statistical ensemble of particle- oscillator systems. The initial oscillator states chosen 
were j = 0, 1, and 2. As seen from appendix E, most surface atoms envisioned as oscil- 
lators at surface temperatures in the range from 5 to 300 K would be in these low states. 

Table I gives the normalized, real probability current densities and the three density 
sums, for a typical set of input parameters shown in figure 3, for the helium- tungsten 
combination. 


15 



oo 



Figure 3. - The input parameters for a sample case of helium-tungsten combina- 
tion. Mass ratio, y, 0.2087; initial oscillator state, j, 1; fundamental frequency 
of oscillator, oj , 1. 31x 10*3 per second. Current densities presented in table I. 


TABLE I. - CURRENT DENSITIES AND DENSITY SUMS FOR SAMPLE CASE OF 
HELIUM- TUNGSTEN COMBINATION 

[input parameters shown in fig. 3. Dimensionless energy, £q, 0. 50; dimension- 
less well width, n, 20.76; dimensionless well depth, 0, 0.7045; mass ratio, y, 

0. 2087; initial oscillator state, j, 1; fundamental frequency of oscillator, co, 

1. 3lxl0 13 sec -1 . C-state sum = D-state sum = 1. 5519; B-state sum 
= 1 . 0000 .] 


Final 

Well current density 

Transition probability, 

cr/v 

oscillator 

energy 

s \ 1/2 z 

A 0 + j - f + /A 1/2 | 2 

level, 

f 

l «. ) 1,1 

l «. J 1,1 


0 

0.0095 

0. 1777 

0. 1831 

1 

1. 3489 

1. 1787 

.8169 

2 

. 1935 

. 1955 

Imaginary 

3 

Imaginary 

Imaginary 

Imaginary 


If the initial incoming particle energy e Q is allowed to vary, the C- and D-state 
sums vary with Cq, as shown in figures 4(a), (b), and (c) for j = 0, 1, and 2, respec- 
tively. Each figure gives results for the two well widths mentioned previously. As 
shown in appendix C, these sums need not be equal to 1, and a possible physical explan- 
ation for their variance from 1 is now offered. 

Requiring the C- and D-state sums to be equal to 1 would apparently be requiring, in 
a physical sense, that the same number of particles be found moving to the right and to 
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Normalized probability current density sum in well 




(b) Initial oscillator state, j, 1. 



Figure 4. - Probability current density sum in well for helium-tungsten combination. Fundamental frequency 
of oscillator, a), 1. 3 lx 10 13 per second; mass ratio, y, 0.2087; dimensionless well depth, p, 0.7045. 
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the left in the well region as are coming in toward the well from + °° with initial momen- 
tum p(j). However, both reflection and transmission can occur at the barrier r - R = b, 
regardless of particle energy or direction of impingement. Also, definite reflection is 
occurring at the other end of the well, the infinite wall. For certain choices of input pa- 
rameters, all particles coming toward the well would not be expected to enter the well. 

In these cases, the C- and D-state sums, envisioned as "mass counters, ” would not 
"count" these particles at all and might be expected to be less than 1. Now, those par- 
ticles that do enter the well region can suffer multiple collisions and accompanying 
changes of their energy. One interpretation is that these well sums may "count” some 
particles more than once. The same particle could conceivably contribute to current den- 
sities for various component states in ^(r, R), associated with particle motion in the 
same direction (left or right), as the result of multiple collisions plus accompanying en- 
ergy transfers. Therefore, the sums might be expected to be greater than 1. Only the 
complete mathematical solution can reveal which situation predominates for a given set 
of input parameters and whether the sums will be less than, equal to, or greater than 1. 

However, for low incident particle energies, initial reflection from r - R = b might 
be expected to predominate, as it does for a low energy particle impinging on a similar 
barrier fixed in space. Then, the well sums would be less than 1, and this is indicated 
in figure 4. Also, for high incident particle energies (eg » / 3 , the well depth), the well 
sums might be expected to approach the value they would have if the well were not there; 
that is, 1, and this is indicated in figure 4. 

Figure 4 shows that the current density sum in the well has a maximum value in the 
€q interval 0. 3 < eg < 0. 5. If this sum is a measure of the frequency of multiple colli- 
sions within the well, a helium atom with eg in this interval has the best chance of being 
trapped, in the "apparent" sense of trapping. 

Results show that neon-tungsten well sums have more local maximums over the £q 

range of figure 4 than helium-tungsten well sums. And the magnitudes of these maximums 

are much greater for neon than for helium. For example, for the neon- tungsten combina- 

13 

tion, the following input parameters were chosen: y = 0. 4686, u> = 1. 31x10 per second, 

/3 = 3. 17025, p. = 20. 76, and j = 0. Again, the energy eg was varied from 0. 025 to 
1. 90, and well dimensions were based on Goodman's Morse potential data (ref. 11). The 
input for the corresponding helium- tungsten case was the same as for neon-tungsten, ex- 
cept with j3 = 0. 7045 and y = 0. 2087. 

The well density sum for neon- tungsten had about four local maximums in this €q 
range (figure not shown). The values of these maximums varied from 3. 1 to 9. 2. The 
corresponding helium- tungsten case is presented in figure 4(a). This figure shows that 
the well sum for ju = 20. 76 has only one local maximum, which occurs at €q - 0. 45. 

Its value is about 1. 65. The greater number of local maximums over a given eg range 
and greater magnitudes of these local maximums for neon may suggest higher sticking 
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probabilities for neon than for helium on tungsten. Goodman also predicted this differ- 
ence in sticking probabilities for helium and neon by using classical mechanics (ref. 11). 
Unfortunately, to our knowledge, experimental trapping data for monoenergetic helium 
and neon beams incident on tungsten does not exist. 


Relative Magnitudes of Cf-State arid Df-State Densities 

An interesting and, as yet, unexplained result occurs whenever the input parameters 
(for helium or neon on tungsten) are such that there exists only one possible stationary 
Bj state. This state is the same as the initial incoming state with the opposite sign for 
particle velocity; that is, the particle can only leave the well with an energy equal to its 
initial incoming energy. From energy considerations, this can happen only when j = 0 
and Cq < 1. When this occurs, not only are the C- and D-state sums equal, but also the 
individual corresponding C^- and D^- state densities are equal (to nearly eight- significant- 
place accuracy of the computer printout); that is, 

M lc , 2= M| Df |2 

m m 1 

for all values of y, ju, and j3 examined thus far. When eg is varied to slightly above 1, 
and then two Bj states are possible (Bq and Bj), this term-by-term equality is de- 
stroyed, but the C- and D-state sums are still equal. 

The equality is destroyed in a somewhat predictable way, though. In the statistical 
interpretation, because of the direction of particle motion associated with and 
states, particles in states could include some that have not "struck'’ the infinite wall, 
whereas all particles in Dj states have "struck" the infinite wall at least once. If most 
of the energy transfer occurs through collisions with the infinite wall, the set of C^- state 
densities might be expected to indicate a higher probability of the initial state j than the 
Dj- state densities. Indeed, the density is always larger than the density, except 
in cases where the term-by-term equality mentioned previously holds (see e. g. , table I). 


Bound- State Sums 

The bound-state sums defined previously were not extensively investigated. The 
physical phenomenon of sticking implies the presence of the gas particle near the surface, 
and it was believed that all states in the well (and their associated current densities) 
should be of interest - not just those corresponding to particle energies less than zero. 
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A particle in the ell can be reflected at the barrier r - R = b back toward the wall re- 
gardless of its total energy. Available results for the helium- tungsten combination indi- 
cate that these two bound-state sums are not, in general, equal, but that they increase 
with increasing well depth jS, as might be expected. 

Energy Accommodation Coefficient Results 

The comparison between <* calc and experimental data for « ac , for helium on tung- 
sten, is shown in figure 5. The experimental curve was drawn from data compiled from 



Figure 5. - Energy accommodation coefficient of helium on clean, 
cold tungsten. Experimental curve taken from data compiled by 
Trilling (ref. 16). Gas temperature assumed equal to wall temper- 
ature (see discussion following eq. (22)). Mass ratio, y, 0. 2087; 
fundamental frequency of oscillator, u), 1. 31xl0 13 per second; 
dimensionless well width, \i, 20. 76; dimensionless well depth, p, 

0. 7045. 


several sources (Silvernail, Thomas, Tho-Nhan, Wachman) and presented by Trilling 
(ref. 16). To enable a comparison between ® ca j cJ for a monoenergetic gas beam and 
a exp f° r a Maxwell -Boltzmann gas, the approximate correlation ^ kT — €qHw was again 
made. (See the section Results for Helium and Neon on Tungsten). Note that none of the 
input parameters (in particular (3 and n) were specially chosen to give the best agree- 


ment between a 


and a , 


calc exp' 

For the temperature range from 20 to 120 K, the simple theory underestimates the 
accommodation coefficient at the lower temperatures and overestimates it at the higher 
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temperatures. In particular, the simple theory yields a limit of zero for a ca ± c , as the 
surface and gas temperatures approach zero. As the limit of zero temperature is ap- 
proached, practically all the oscillators are in the ground state (3 = 0). The gas particle 
approaches the oscillator with energy £q « hco, and no phonons can be exchanged after 
the gas particle returns to +°° after interacting with an oscillator in state 3 = 0 . There- 
fore, E in = E qu ^ in equation (23), and ol = o!q = 0. Therefore, c7 calc approaches zero. 
Goodman's classical theory (ref. 11) has a low temperature limit of 1, and a minimum 
with varying temperature for helium- tungsten which is supported by recent experiments 
(fig. 5). 

Goodman's theory yields a high temperature limit for the effective accommodation 
coefficient of a single gas atom of 4(m/M)/[l + (m/M)]^, which equals 0. 0833 for the 
helium-tungsten combination. This value seems to be exceeded by a ca i c at the higher 
temperatures of figure 5. However, it must be noted that this familiar classical hard 
sphere limit of 4(m/M)/[l + (m/M)] is a limit only as the incident gas particle energy 
(temperature) approaches +°°. The energy (temperature) of the solid atom is still taken 
as zero in this limit. On the other hand, the temperature of figure 5 is (approximately) 
both the surface and incident gas temperature, in accordance with the limit definition in 
equation (22). The high temperature limit of a ca i c in figure 5 has not been ascertained 
because the computer cannot calculate accurate transition probabilities (and hence E qu ^. 
in eq. (23)) when Cq and 3 become large. 

There are several possible reasons for the discrepancy in figure 5. Calculations 
were performed only for one fundamental frequency, where actually, a frequency spec- 
trum (e. g. , Debye) should be considered. The calculation of ct^ could be repeated for 
many frequencies, and by weighting the results in accordance with the relative occurrence 
of those frequencies in the Debye spectrum, possibly a more accurate accommodation 
coefficient could be derived. Also, the calculation could be repeated for many incoming 
particle energies, and the results weighted in accordance with the relative occurrence of 
those energies in a one-dimensional Maxwell- Boltzmann energy distribution, to give a 
theoretically more accurate gas-temperature - gas- energy correlation. Such calculations 
have not as yet been performed because they require excessive computer time. 

Furthermore, the trouble could be inherent in the model. Allowing the gas particles 
to make only direct, head-on, one- dimensional collisions with single surface layer parti- 
cles, with no energy transfer to and from the neighboring lattice particles, could quite 
possibly yield misleading values for the accommodation coefficient. Also, the assumed 
interaction potential is a rather poor approximation to the more realistic Lennard- Jones 
or Morse potentials. Attempts to remedy most of these model deficiencies lead to a con- 
siderable increase in the complexity of a quantum- mechanical solution. 


Since the calculation of <* calc 
totic region, 


involves only the transition probabilities in the asymp- 
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r e n + 5 - f 


1/2 




the importance of including the well in the model for a study of accommodation coeffi- 
cients may be questioned. That the well presence significantly affects such probabilities, 
and hence ® ca i c > is shown in table n. In each of the two sample cases (eg = 0. 400 and 
€q = 1. 500), the probability that the oscillator is found in its initial state (j = f = 1) when 
the particle is approaching +<*> after the interaction is definitely smaller when the well is 
present than with no well at all. 


TABLE n. - EFFECT OF WELL PRESENCE ON TRANSITION 
PROBABILITIES FOR HELIUM ON TUNGSTEN 


[Mass ratio, y , 0. 2087; dimensionless well width, p, 20. 76; 
initial oscillator state, j, 1.] 


Final 

Dimensionless 

Initial 

Transition probability, 

oscillator 

energy 

well depth, 

P 

energy of 
particle, 


level, 



0 

V ‘0 / ' 

f 




0 

0 

0.400 

0.0628 

1 

0 



.9372 

2 

0 



Imaginary 

0 

.07045 



0. 1703 

1 

.07045 



.8297 

2 

.07045 



Imaginary 

0 

0 

1. 500 

0. 1502 

i 

0 



.7156 

2 

0 



. 1342 

3 

0 



Imaginary 

0 

. 07045 



0. 1649 

1 

. 07045 



. 5254 

2 

.07045 



.3097 

3 

.07045 



Imaginary 
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SUMMARY OF RESULTS 


The Schrodinger equation was solved numerically for a one-dimensional collision be- 
tween a harmonic oscillator and a particle with a rectangular trapping potential. General 
solutions were sums of states, each a product of a particle wave function and an oscillator 
wave function. 

A correlation was suggested between probability current densities for states within 
the well and sticking of gas atoms on surfaces, where sticking is regarded as arising from 
multiple collisions. Energy accommodation coefficients were formulated and compared 
with helium- tungsten experiments for the temperature range 20 to 120 K. The theory un- 
derestimates the accommodation coefficient at lower temperatures and overestimates it 
at higher temperatures. Results indicate that only approximate qualitative predictions 
of the coefficient are possible with this simplified model. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 11, 1968, 

124-09-19-01-22. 
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APPENDIX A 


SYMBOLS 


Aj, Bj, Cj, Dj coefficients in ^(r, R) and 
* 2 ( r ’ R ) 

a constant in Morse potential, 

X-l. „-l 


m 


H n («R) 


b 

D 

E. 


m 


E 


osc 


J out 


Er 


E(1),E(2) 


F, K 


well width, m 

constant in Morse potential, 
kcal/mole; J/mole 

initial particle energy, 
e 0 Ko>, J 

harmonic oscillator ener- 
gy, (n + n = 0, 1, 

average outcoming particle 
energy, calculated in 
text and appendix E, J 

average energy of inci- 
dent particle in defini- 
tion of Of , J 
ac’ 

average surface molecular 

energy, in definition of 

a , J 
ac’ 

average energy of gas par- 
ticles leaving surface, in 
definition of a , J 

SIC 

Schrodinger equation ener- 
gies for regions 1 and 2, 
respectively, J 

cutoff indices 

largest value of f for 
which jp-^f)] is real 


“2 


K 

k 

M 

m 

P(J) 


P M 


largest value of f for 
which is real 

Her mite polynomial 

h/27r, where h = Planck's 

_ Q4 

constant = 6. 6256x10 
J-sec 

derived quantity, ye q^ 2 

derived quantity, 
y(e 0 + j - f) 1/2 

derived quantity, 
y(e Q + j - f + /3) 1 / 2 

positive integers denoting 
initial and final oscilla- 
tor states, respectively, 

0 1 00 
V y x j f • • 

summation index 


Boltzmann constant, 

1. 38054X10' 23 J/K 

mass of oscillator, kg 

mass of particle, kg 


fraction of oscillators in 
state j, from Boltzmann 


exp 


/ m<»\ 

l-exp/- Rw \l 

V kT s/ 

Hi] 


oscillator transition prob- 
ability, 


+ j -f 


a/2 


B. 


'fl 
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P(j) 

Pj(f) 



expansion coefficient, defined 
in text, l = 0, 1, 2 

initial particle momentum, 
(kg)(m)/sec 

particle momentum, in region 
1 when oscillator is in state 
f, (kg)(m)/sec 

denotes choice of positive root, 
real or imaginary, when 
Pj(f) is calculated from 
eq. (5) 


P2< f ) 



particle momentum, in region 
2, when oscillator is in 
state f, (kg)(m)/sec 

denotes choice of positive 
root, real or imaginary, 
when ? 2 (f) is calculated 
from eq. (6) 


Q K (hj , f) expansion coefficient, defined 
in text, l = 0, 1, 2 


R 


r' 

V 


coordinate of oscillator, A; 
m 


reduced particle coordinate, 
C 0 


r' - X n + * n - , A; m 


coordinate of particle, A; m 

gradient operator with respect 
to r 

r Q nominal radius of oscillator, 

m 

Tp nominal radius of particle, m 

S total current density 

SrOij , f) expansion coefficient, defined 
in appendix C, l = 0, 1, 2 


S l- S 2 

T 

W 1 * 


V 


Mm 


v(r* - R) 


X 

X, 


a 


a 


ac 


total particle current densities 
in regions 1 and 2, respec- 
tively 

temperature, K 

expansion coefficient, defined 
in appendix C, l = 0, 1, 2 

surface temperature, K 

temperature of incident parti- 
cles, in definition of a , K 

temperature of surface mole- 
cules in definition of a , K 

temperature of gas particles 
leaving surface, in definition 

of “ac’ K 

rectangular interaction poten- 
tial, J 

well depth in V Mm , J 

Morse potential, 

Djexp£-2a(r* - R - Xq)J - 

2 exp^a(r’ - R - X 0 )Jj, 
kcal/mole; J/mole 

dimensionless coordinate of 
oscillator, aR 

equilibrium separation in 
Morse potential, A; m 

dimensionless reduced parti- 
cle coordinate, ar 

reciprocal of characteristic 
harmonic oscillator length, 




m 


-1 


energy accommodation coeffi- 
/E„ - E, 

cient, lim 


E,-E 0 , 


E 0 - E 1 1 
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a calc thermal average of cr values, 

2 p Q) a i 

3=0 ] 

Q! ___ experimental value of a for 

6Xp 3.C, 

helium on tungsten, taken 
from data compiled by 
Trilling (ref. 16) 

a!j microscopic energy accommo- 

dation coefficient, 

(E. - E .)/(E. - E ) 

' in out' /v in osc' 

/3 dimensionless well depth, 

Vj/hw 

y mass ratio, ( 2 m/M)^ 

6 ^ Kronecker delta: 1, if f = j; 

0 , if f i j 

€ q dimensionless energy, 

[p(j)] 2 / 2 mfio> 


e j dimensionless energy, 

jp 1 (f)J 2 / 2 mKa) 

eg dimensionless energy, 

[p 2 ( f )] 2 / 2 mKu) 

#0 Debye temperature, K 

H dimensionless well width, ab 

£ normalization constant in 

* n <r), [a/fr^V)] 1 / 2 

^(r, R), Schrodinger equation solutions 
SI'gCr, R) for regions 1 and 2 , respec- 
tively 

^ n (R) harmonic oscillator wave func- 

tion, £ n exp (-j a 2 R 2 jH n (aR); 
n = 0 , 1 , ... oo 

co fundamental frequency of oscil- 

lator, sec - * 
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APPENDIX B 


RELATION BETWEEN MORSE AND RECTANGULAR WELL POTENTIALS 

Equilibrium occurs in the Morse potential when the coordinates of the oscillator and 
particle are separated by the distance X^; that is,, when r* - R = Xq. Then, the Morse 
potential (fig. 6) can be written as 

v(r' - R) = D|exp[-2a(r’ - R - X Q )] - 2 exp[-a(r’ - R - X Q )]| (Bl) 

The parameter D is the depth of the Morse potential well, and the reciprocal of the pa- 
rameter a is a measure of the width of the well. 

A rectangular well, with one side as an infinite wall, will be used to approximate the 
Morse potential. During development of this approximation, for convenience in subse- 
quent calculations, a new coordinate r, linearly related to r', will be introduced. 

The coordinate of the infinite potential wall on the (r' - R) scale of figure 6 will be 
fixed arbitrarily at 

< r ’ - R >wall - *0 ~ < B2 > 

If a new coordinate is defined 



separation, X Q 

Figure 6. -Morse potential v(r' - R) and rectangular well potential V Mm . 
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equation (B2) can be written as 


r’ - R - X Q + 


In 2 


- ( r - R ) wa ll “ 0 


'wall 


(B3) 


The other end of the rectangular well, the potential step, will be fixed at 

< r '- R >step= X 0-lT + b (B4) 

a 

where b is the rectangular well width (see fig. 6). Again, using the coordinate r gives 
equation (B4) in the form 

(*■' - R - x 0 + — ) =< r - R > s tep = b < B5 > 

\ a /step 

The transformation to coordinates r and R is not absolutely necessary. The alge- 
braic form of some subsequent equations is slightly simpler if one end of the rectangular 
well occurs at a point where the difference of the coordinates chosen is zero 
(r - R) wa u = 0 . Thus, the interaction potential V Mm becomes that shown in figure 2. 

For this initial study, one value of well depth Vj = D was chosen, and the well width 
b was treated as a variable parameter. For the helium- tungsten combination, two values 
of b were arbitrarily chosen, b = (2 In 2)/a and b = (3. 5 In 2)/a. Using^Good man's val- 
ues of a = 1. 3 per angstrom and D = 0. 14 kilocalorie per mole (5. 86x10 J/mole) for 
this combination (ref. 11), o> = 1. 31xl0 13 per second (using appendix B of ref. 11), and 
equation (11), gives the well depth /3 as 0. 7045 and the two well widths as ju = 20. 76 
and 36. 33, respectively. 
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APPENDIX C 


REDUCTION OF BOUNDARY-CONDITION EQUATIONS 


Consider equation (8). The normalized oscillator wave functions are given by 


VX) = «„exp(-i X 2 ) H n (X) 


X = a R 


(Cl) 


where 


£n = 


a 


1 1/2 


11 \7T 1/2 2 n n:> 


and 




dX = a 


Using the relation (d/dX)H n (X) = 2nH n _j(X), one can show that 

® = [ffi ^" (X> ] ^ = “ «n + < C2 > 

Substituting equations (3) and (4) into equation (8) and using equations (11) and (C2) yield 
the following equation: 


L. exp[-ih 0 (X + M)]«|j exp^y.XHjfX) + 2JH j _ 1 (X)] 


+ ^ B f exp 
f=0 


/ 2 \ 00 

1+ihjpC + M)]«4 f fxH f (X) + 2fH f _ 1 (X)Jexp|- — J = £ [c f exp[-ih 2 (X + ju)] 

\ 2 / f=0 L 


+ D f expJ+ih 2 (X + ji)]|a^ f exp ^ j|xH f (X) + 2fH f l (X)j 


(C3) 


Now consider the combination of the form ijjjHj_j(X) in equation (C3). From the defini- 


tion of £ n , it is apparent that 
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- 11/2 




a 


j-1 


or 


L^-VdJ 


^ j = (2jr 1/2 ^_ 1 


= (2j) l/2 ^ 


3 


(C4a) 


(C4b) 


Therefore, 


/ 2\ 

A. exp|jih 0 (X + pj]o' exp f-^-j (2j)" 1 / 2 (2j)^._ 1 H._ 1 pC) 


= A. expJj-ih Q (X + /i)]«(2j) 1,/2 i//._ ^X) (C5) 

Us ing equations (C4) and (C5), with the possible substitution of f for j in equa- 
tions (C4), reduces equation (C3) to 

Aj exp£-ihQ(X + M)][-« x ^ (x) + a(2j) l/2 ^_ 1 (X)j 

+ Yj B f expJ+ih^X + m)] |jaXi// f (X) + (2f) 1,/2 Q'i// f _ 1 ( x )J = f^|c f expf-il^X + /i)] 


+ D f exp[+ih 2 (X + M )]|^-ffX^ f (X) + (2f) l/2 a l /, f _ 1 (X)J 


Some terms in equation (C6) are expanded as follows: 


(C6) 


exp(+ih z X)X i// f (X) = £ T K (h z ,f)^ K (X) 
K=0 


where l = 0, 1, or 2, and 


a 


/oo 

^ exp(ih z X)X^/ f (X)^ K (X)dX 

= (f)i;[" QK(hi ’ f) ] 


(C7) 


(C 8) 
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Using equations (C8) and (15) shows, after some rearrangement, that 

*m 

^2 1 , . I /u \\ \ / i\to^/u \-2t 


TK<h ' ,f)= (^l) exp 


_/V\ h f+K/ i f+K+l\ (\ _ f + k\ \ X (- 1) 2 (hp 
\2/ 1 ' ' \2 h z // y t:(K - t)I(f - t): 

- J L t=o 


m 


■E 

t=0 


(- l) t+1 2 t+1 t(hj)" 2t-1 
t:(K - t):(f - 1): 


(C9) 


where 1=0, 1, or 2, and t is the smaller of the values f or K. Some terms in 
equation (C6) are also expanded as follows: 


exp(-ih z X)X 1 // f (X) = £ S K (h z ,f)^ K (X) (CIO) 

K=0 

with Sg-(h^, f) having the same algebraic form as T^(h^, f), with +i replaced by -i. 

Now, using equations (C7), (CIO), (13), and (14) gives equation (C6), after inter- 
changing orders of summation and letting 


6 . =fl ^ f = 1 

fj \o if f^j 


in the form: 

r 


E %( x )i 

K=0 


2 A i 6 fj|exp(-ih 0 M)S K (h 0 , j) + (2j) 1//2 P K (h 0 , j - l)exp(-ih 0 p)j 
f=0 L J 

QO 

+ 2 B f exp(ihip)j^T K (h 1 ,f) + (2f) 1//2 Q K (h 1; f --1)J 

OQ 

+ V C f exp(-ih 2 p)[s K (h 2 ,f) - (2f) 1/2 P K (h 2 ,f - 1)] 


£ D f exp(ih 2 p)[T K (h 2 ,f) - (2f) 1/2 Q K (h 2 ,f - 1) 


>=0 (Cll) 


j 
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Now, once again, the completeness of the set |^ K (X)| immediately yields equation (18) 
from equation (Cll). A similar procedure leads from equation (9) to equation (19), and 
from equation (10) to equation (20); no new expansion coefficients are involved. 

Because of the presence of h^ in some denominators in the algebraic forms of the 
expansion coefficients, the computer cannot evaluate the coefficients in these forms when 
hj = 0. It is easiest then to return to the original definitions of the coefficients and to use 
the following equations, when h^ = 0 (ref. 17): 


T k ( 0, f) = S K (0, f) = 1 ^ i//*(X)X^ K (X)dX 


= if K = f + 1 


/f\i / 2 

-Q UK - t - 1 


(C12) 


= 0 Otherwise 


Q k (°, f) = P K (0, f) =- f **(X)* K (X)dX 

Of J- oo A 


= 1 if K - f 


Otherwise 


(C 13) 
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APPENDIX D 


TOTAL CURRENT DENSITY CALCULATIONS FOR REGIONS 1 AND 2 


Using ^(r, R) as given in equation (3) shows that 



Note that'[pj(f)J is now either real or pure imaginary, so that its complex conjugate 
is not, in general, equal to itself. Using equation (Dl) and the orthonormality of the set 



If r is taken large enough, the last exponential in equation (D2) 



will be equal to 1 if [Pj(f)] is r eal and will approach zero if p^(f) = i|p^(f)|. Doing this 
and noting that the third term in equation (D2) is simply the complex conjugate of the sec- 
ond term result in 
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s 


1 



(D3) 


where is the greatest value of f for which [p 1 (f )] is real. Equation (D3) means that 
the B-state sum must equal 1. 

Using ^ 2 (r, R ) as given in equation (4) and the orthonormality of the set |V f (R)j in 
equation (29) yields 


S 2 = 0 = — Im 
m 


/ 7 l 1 B j K 

f=0 V ^ ^ ^ 


l [f2 (f 3 ] 


+ D*C f £M 


K 

v J 


>exp< 


+ c*d J 


i[p 2 ( f jfl Ji(p 2 ^ f ) 


K 

v. J 


+ |D f | 2 ^ 


P[P 2 (f 3 


>exp< 



(D4) 


Next, the infinite sum is broken up and ^ is introduced as the greatest value of f for 
which jp 2 (f)] is real: 


S 0 = 0 = — Im 
1 m 


*2 

Z 

Lf=0 
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Equation (D5) implies that the C- and D-state sums must be equal, but nothing is said 
in equation (D5) about the magnitude of their common value. This value is a complicated 
function of the input parameters A^, £q, j, y, (3, and p. 
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APPENDIX E 


SAMPLE CALCULATION FOR ENERGY ACCOMMODATION COEFFICIENT 

The helium- tungsten combination is considered once again with y = 0. 2087, 
co = 1. 31X10 13 per second, p. = 20. 76, /3 = 0. 7045, and 6 q = 0. 225, corresponding to a 
gas temperature of approximately 45 K. This temperature was the apparent crossover 
point of theoretical and experimental curves for ct (see fig. 5). The surface temper- 
ature T will be assumed to be 45 K also. A difference of a few degrees would only 
slightly repopulate the surface layer oscillators. 

For j = 0: 


P 0~0 = L 000 

All other transitions are energetically impossible in the asymptotic region. Therefore, 

Ejn = 0. 225 Ka> 


E osc = (° + |f" 


E out = (1.0000X0. 225)Ku> 


E. - E . 

out = 0 „ 

E iu- E osc 


For j = 1: 


P l-0 = °* 0930 
P l-1 = °* 9070 


Dropping the hco yields 


E in = 0. 225 

E = 1 + i = 1. 500 
osc 2 
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E out = (°* 9070)(0. 225) + (0. 0930)(1. 225) 

« 0. 318 

Note that if the oscillator went from j = 1 to f = 0, the particle energy must have 
changed from 0. 225 to 1. 225, so that' 1. 225 is weighted by P-j^q in E^: 

a _ 0.225 - 0.318 
1 0. 225 - 1. 5 

« 0.0730 


For j = 2: 


E out 


P 2 — 0 = °* 0110 

P 2-l = °* 1900 
P 2 _ 2 = 0.7991 

E. n = 0. 225 
E osc = 2 * 5 

= (0. 7991)(0. 225) + (0. 1900)(1. 225) + (0. 0110)(2. 225) 
0.437 


Therefore, 


For j = 3: 


= 0. 0931 

P Q = 0. 0010 

P 3-l = °* 0335 
P 3_2 = °* 2843 
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= 0. 6812 


P 3-3 
E. n = 0. 225 

E = 3. 5 
, osc 

E out = (°* 6812)(0. 225) + (0. 2843)(1. 225) + (0. 0335)(2. 225) + (0. 0010)(3. 225) 
« 0. 5791 


Therefore, 


Ofg = 0. 1080 

Now, with T g a* 45 K and 6 = Kco/k =* 100 K, 



= (0. 892)(0. 1080) j 


P(0) e* 0. 892 
P(l) « 0. 0962 
P(2) « 0.0104 
P(3) « 0.001122 

Therefore, 

“calc “S Pa) “i“ 0, 00812 


The relative smallness of P(j)a!., for j = 3 and higher, makes possible summing to 
about j = 3 in the expression for ® ca i c . This holds true for the entire temperature 
range of figure 5. 
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